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(54) Body Motion Compensation 
Filter 

(57) A filter for use In determining the 
vertical displacement of a body in 
water and more particularly for 
surveying underwater terrain from a 
towed sensor comprises a first portion 
receiving, 52, a position signal from a 
pressure transducer and having a 
transfer function which is substantially 
the product of a predetermined 
constant and a first real rational 
polynomial function P v (s)/Q(s); and a 
second portion receiving at 54, a 
vertical acceleration signal from an 
accelerometer and having a transfer 



function which is substantially the 
product of a predetermined constant 
and a second real rational polynomial 
function SP 2 (s)/Q(s). The output 
signals from the two portions are 
summed at 70. The degree of Q(s) is 
an integer n greater than 2, and that of 
P t (s) and of P 2 (s) is less than n. The 
magnitude of P t (s)/Q{s) is to be less 
than 1 for s=jw where w is any radial 
frequency within the pressure 
transducer noise spectrum. While the 
zeroes of Q(s) are such that 1/Q(s) has 
a 20n db/decade roll-off in the said 
spectrum P,(s) and P 2 (s) are 
constrained to satisfy the relationship 
P 1 (s)+s 3 P 2 <s)=Q(s). 
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SPECIFICATION 

Body Motion Compensation Filter 

The invention relates to an electronic filter for use in determining the position of a body travelling 
through water relative to a datum such as mean sea level. 
5 The invention has specific application to an underwater seismic system for use in surveying the 5 
terrain under a body of water. Such systems commonly include a towed body, generally referred to as a 
"fish", which is provided with a sound source and hydrophones or like sound sensors. Very briefly, a 
sound pulse is generated and directed towards the bottom terrain and a reflected pulse Is detected by 
the hydrophones. From the delay between generation of the sound pulse and detection of the reflected 

1 0 pulse, an estimate of the distance to the reflecting layer of the terrain below may be obtained. By 10 
towing the fish at a relatively constant speed and by sampling the distance to the bottom terrain at 
regular intervals, it is possible with a chart recorder to produce a profile of the bottom terrain. 

The quality of the profile obtained may be seriously degraded by vertical movement of the fish 
during the sampling process. Upward movement of the fish makes the bottom terrain appear farther 

1 5 from the datum; downward movement makes the terrain appear closer. A very significant factor in this 1 5 
respect is heaving of the fish as the towing vessel is buffeted by surface waves. This heaving motion 
translates onto the chart recorder as an oscillatory pattern which tends to reduce the resolution and 
usefulness of the profile. U.S. Patent No. 4,01 9,356 to Hutchins discloses apparatus for compensating 
for such heaving motion, including pressure and acceleration transducers for use in determining the 

20 instantaneous vertical displacement of the fish with respect to the mean sea level. 20 
Pressure and acceleration transducers presently used are subject to noise problems. This noise 
tends to degrade the quality of chart recordings in a manner not unlike that of the heaving motion 
experienced by the fish. 

More specifically, the pressure transducer used will generally derive from pressure measurements 

25 a position signal indicative of the vertical! displacement of the fish with respect to a datum such as 25 
mean sea level. The positive signal is likely to be subject to noise owing to variations of pressure 
caused by surface wave action. Since this same wave action buffets the towing vessel!, the pressure 
variations tends generally to have a frequency range similar to that of the heaving motion. This overlap 
in spectra makes it difficult to attenuate pressure transducer noise without affecting the portion of the 

30 position signal reflecting the heaving motion of the fish. 30 
The accelerometer used will generally produce an acceleration signal proportional to the 
acceleration of the fish along a specific axis through the accelerometer. The axis will normally be 
arranged to be in a substantially vertical position when the fish is moving smoothly through water so 
that the acceleration signal is substantially indicative of the vertical acceleration of the fish. Such an 

35 accelerometer can generally be expected to have a zero frequency offset corresponding to the effect of 35 
gravity on the accelerometer. Because the accelerometer axis tends to tilt from a vertical position as 
the accelerometer moves with the pitching and rolling of the fish, the gravity offset will at any instant 
be attenuated by the cosine of the angle by which the axis is tilted from a vertical position: the gravity 
offset is in effect modulated by the pitching and rolling of the fish. 

40 The modulated gravity offset is an accelerometer noise signal which can be expected to have a 40 
very broad frequency range and a significant zero frequency component owing to the fact that over the 
typical range of angles formed by the accelerometer axis with respect to a vertical position the cosine 
function provides a positive result. The zero frequency noise component can in practice be the principal 
portion of the accelerometer noise and can be different to remove without distortion of the vertical 

45 acceleration signal and ultimately of the displacement signal to be produced. 45 

The frequency spectra of the various signals will vary depending on the physical circumstances in 
which the fish and the towing vessel are located; however, it has been found practical for design 
purposes to assume that the period of the heaving motion will be in the order of 6 to 1 0 seconds with a 
resultant frequency spectrum centered about . 1 25 Hz. As mentioned above, the frequency spectrum of 

50 the pressure transducer noise can tend to overlap of the heaving motion. 50 
A filter constructed according to the invention effectively comprises first and second filter 
portions those output signals are added by an electronic summer. The first filter portion is intended to 
receive and filter the position signal produced by the pressure transducer; the second filter portion, the 
acceleration signal produced by the accelerometer. The transfer functions of the first and second filter 

55 portions are chosen to attenuate the noise characteristic of the signals respectively received, and are 55 
related in a manner described more fully below that permits a signal indicative of the body's vertical 
displacement to be obtained by summing the first and second filter portion output signals. 

The words "transfer function" as used in this disclosure and in the appended claims denote the 
Laplace transform of the Impuse response of a filter or like electronic circuit. The transfer functions 
60 dealt with will be real, rational, polynomial functions in the complex Laplace variable s. Transfer 60 
function "zeroes" and "poles" will consequently be the values of the Laplace variable s at which the 
numerator and denominator polynomials of the transfer function respectively assume values of 0. It will 
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According to the invention, the first and second filter portions will respectively have transfer 
functions T,(s) and T 2 (s), of the following general form: 

M P,(s) 

TVs)- 



*1 Q(s) 
1 sP 2 (s) 



T 2 (s)=- 

Ka Q(s) 

5 In the above, M is a gain factor which may be freely chosen; K, is the accelerometer gain, the constant 5 
of proportionality relating the body's acceleration to the accelerometer output signal; and, K 2 is the 
pressure transducer gain, the constant of proportionality relating the vertical displacement of the body 
to the pressure transducer signal. P-{s), P 2 (s) and Q(s) are polynomial functions in the complex variable 
s. Because the transfer functions T t (s} and T 2 (s) are constrained to be real, rational!, polynomial 
1 0 functions, the coefficients of P^s), P 2 (s) and Q(s) will of necessity be real numbers. The degree of Q{s) \ q 
is an integer n greater than 2, and the degree of each of P,(s) and P 2 (s) is less than n. P/s) and P-(s) are 
further constrained to satisfy the relationship: 

P 1 (s)+s 3 P 2 (s)=Q{s) 

If the Laplace transform of the body's vertical displacement with respect to the datum is a 
1 5 function Z(s) and consequently that of the acceleration of the body's s*Z(s), then the Laplace transform -a 5 
of the filter output signal E(s) will be given by: %# 

Simple algebraic manipulation of the expression above leads to the expression: 

E{s)=MZ{s) 

20 Thus, a signal proportional to the vertical displacement of the fish will appear undistorted at the filter 20 
output. 

The pressure transducer noise received by the first filter portion produces at the output of the 
filter a noise signal whose Laplace transform is the product of the transform of the pressure transducer 
noise and the transfer function T,(s). If the zeroes of Q(s) are chosen so that a transfer function 1/Q(s) 
25 were to have a 20 n db/decade roil-off in the noise frequency range, then, because the degree of P,(s) 25 
is less than the degree of Q(s), the magnitude of the output noise signal corresponding to pressure 
transducer noise may be expected to drop at least 20 db/decade with increasing frequency. If the 
zeroes of P,(s) are now chosen so that the magnitude of P t (s)/Q(s) is less than 1 for s=jw where w is 
any radial frequency in the frequency range of the pressure transducer noise, then an increase in the 
30 signal to pressure transducer noise ratio may be expected. 39 

The forms of P t (s) and Q(s) will be completely specified by the choice of zeroes above. P 2 {s) will 
also be completely specified, owing to the relationship among P.(s), P 2 (s) and Q(s) namely, 
P 1 (s)+s 3 P 2 (s)=Q(s). However, if the zeroes of the polynomial Q(sj are chosen to be non-zero then the 
second filter portion transfer function T 2 (s) will be constrained to have a zero frequency zero. 
35 Consequently, zero frequency components corresponding to the effect of gravity on the accelerometer 35 
will be removed by the second filter portion. 

The invention will be better understood with reference to drawings in which: — 

Fig. 1 is a schematic representation of a filter network which realizes transfer functions related 
according to the invention; 

40 Fig. 2 is a schematic representation of a practical filter constructed according to the invention; 40 
Fig. 3 is a diagrammatic representation of an underwater seismic system in association with 

which a filter constucted according to the invention is intended to be used; and 

Fig. 4 is a schematic representation of the filter of Fig. 2 in use in the underwater seismic system of 

Fig. 3. 

45 Reference is made to Fig. 1 which shows a filter 1 0 having first and second input terminals 1 2, 45 
1 4 for respectively receiving position and acceleration signals as described above, and an output 
terminal 1 6. The filter 1 0 comprises first and second scaling means 1 8, 20; first and second 
differentiation networks 22, 24; a summing network 26; and a filter portion 28 having a transfer 
function 1/Q(s), where Q(s) is a polynomial function of degree n, as described above. 

50 Q(s) will generally have the following form: 50 

Q(s)=a 0 +a 1 s+a 2 s 2 +... +a k s k +... +a„_ 1 s n ~ , +s n 
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by first choosing the zeroes of Q(s) so that the transfer function 1/Q(s) has a 20 n db/decade roll-off in 
the frequency range of the pressure transducer noise, and then producing Q(s) by expanding the 
product of n terms each comprising a different zero subtracted from the complex variable s. With the 
zeroes and coefficients of Q(s) specified in this manner, the transfer function 1/Q(s) may be realized by 
5 using conventional filter design techniques. 5 

The scaling means 1 8, 20 have transfer functions which are respectively a constant M/K % and a 
constant Wl/K 2 , as described above. The purpose of the first and second scaling means 1 8, 20 is to 
compensate for the different scaling provided by the accelerometer and the pressure transducers. The 
transform of the output signal of the first scaling means 1 8 wilt be MZ{s), and that the output signal of 
1 o the second scaling means 20, s 2 MZ(s). The signals will of course be subject to noise. 1 0 

The first differentiation network 22 comprises n-1 unity gain differentiators 30 connected in 
series to produce successive derivatives of the scaled position signal. The second differentiation 
network 24 comprises n-2 unity gaUn differentiators 32 connected in series to produce successive 
derivatives of the scaled acceleration signal; Each of the differentiators 30, 32 has a transfer function 
15 which is simply the complex variable s. r . 15 

The summing network 26 comprises n individual summers 34 which are connected as shown to 
produce a weighted sum of the scaled pressure transducer signal), and the output signals of each of the 
differentiators 30, 32. The weighting used is indicated by constants adjacent to the input terminals of 
the summers 34, the input terminate being generally indicated by arrows contacting the summers 34. 
20 It will be appreciated that the degree n of the polynomial Q(s) has not been specified so that the 20 
filter 1 0 provides very general means for realizing transfer functions (related according to the invention. 
Consequently, the entirety of the filter 10 is not shown and a filter portion 36 is provided in Fig. 1 as 
one of n-5 similar filter portions connected series between points A and points B of the filter 1 0. 

Weighting constants a v and a 2 shown in Fig. 1 are the same a^ a, and a 2 which are 
25 coefficients of the polynomial Q(s). The weighting constants denoted generally by b^ and c b where Ec 25 
assumes values ranging from 3 up to n-1 are related by the relationship b^c^a^ where a h once again 
denotes each of the coefficients of the polynomial Q(s). 

The differentiation network 22 in combination with the summing network 26 constitutes a filter 
portion that realizes a transfer function P^s); the second differentiation network 24 in combination 
30 with the summing network 26 constitutes a filter portion that realizes a transfer function the 30 
polynomial sP 2 (s), as described above. The constants b^ will be chosen to provide the polynomial P,(s) 
with appropriate zero values. The constants c^ will be completely specified by the choie of die 
constants b k . It will now be apparent that the transfer function T^s) is realized between the output 
terminal 1 6 and the input terminal 1 2; and the transfer function T 2 (s), between the output terminal 16 
35 and the input terminal 14. 35 

Although the filter 10 could be realized using frequency limited differentiators for the 
differentiation networks 22, 24 and conventional summing circuits for the summing network 26, in 
practice an equivalent circuit, such as that shown in Fig. 2, is preferably used. Fig. 2 shows a filter 50 
having input terminals 52, 54 for respectively receiving position and acceleration signals as described 
40 above, and an output terminal 56. The filter 50 is intended to realize between the output terminal 56 40 
and the input terminal 52 a transfer function T,{s) as described above, and between output terminal 56 
and input terminal 54 a transfer function T 2 (s), as described above, such that: 

P t (s)=3w 0 s 2 +3wgs+wg; 

P 2 (s)=s; 

45 and 45 

Q(s)=(s+w 0 ) 3 

where w 0 is a positive, real constant. The constant w 0 represents the radial frequency at which a 
transfer function 1/Q(s) begins to have a 60 db/decade roll-off. It will be appreciated that P^s), P 2 (s) 
and Q(s) satisfy the relationship P 1 (s)+s 3 P 2 (s)=Q(s). 

50 The filter 50 comprises four operational amplifiers 58, 60, 62, 64, having inverting terminals 66 50 

and output terminals 68. For the sake of clarity only those impedances connected to the inverting 
terminal 66 and output terminals 68 are shown. The filter 50 also includes first and second summers, 
70, 72 for producing weighted sums of signals, the weighting having magnitude 1 and sign as 
indicated in Fig. 2 beside input terminals marked by arrows contacting the summers 70, 72. 

55 With M, K t , K 2 and Wq, defined as above, and resistor values R 2 and R 3 freely chose, the following 55 

set of component values will enable the filter 50 to realize the transfer functions T,(s) and T 2 (s), as 
defined: 



I 
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MR 3 ! 
R a =3R 2 c i= C 2 = 



K 2 wgR| w 0 R 2 
3W3 



C 3 = 



K 0 w 0 R 2 

The amplifiers 58 and 60 present to the position signal received at the termina! 52 a first low 
pass filter, whose transfer function comprises two poles substantially at s=-w 0 . The amplifier 62 
5 presents to the position signal received at the terminal 52 a first band pass filter, whose transfer 5 
function comprises a zero at s=0 and two poles substantially at s=*-w 0 . The filter 64 effectively 
provides a second low pass filter, whose transfer function comprises a pole substantially at s=— for 
filtering the first band pass filter signal. The second summer 72 produces a signal proportional to the 
weighted sum of the first and second low pass filter signals, and the first band pass filter signal, the 
1 0 weighting being as discussed above. These elements effectively constitute a first filter portion as 10 
described above. 

The amplifier 58 presents to the acceleration signal received at the terminal 54 a second band- 
pass filter, whose transfer function comprises a zero at s=0 and two poles substantially at s=-Wp The 
amplifier 60 effectively constitutes a third low pass filter, whose transfer function comprises a pole at 
11 5 s=— w 0 for filtering the second band-pass filter signal. This effectively constitutes a second filter portion 1 5 
as described above. 

The poles associated with the transfer functions realized by the filter 50 were chosen to be equal, 
negative real constants, namely, -w 0 . It will be appreciated that transfer function poles for filteirs 
constructed according to the invention may generally be chosen from the complex s plane provided 
20 that the requirements specified in the disclosure are observed. 20 
Fig. 3 shows an underwater seismic system 100 comprising a body 102 being towed with a 
cable 104 from a towing vessel 106. A sound pulse 1 08 is shown emanating from the towed body 
102 and travelling towards a reflecting layer 1 10. The rejecting layer 110 causes a reflected sound 
pulse 1 1 2 to be directed upwardly toward the towed body 1 02. Further layers t 14 would also produce 
25 reflected pulses, but the latter are not shown for the sake of clarity. 25 
Fig. 4 shows the filtered 50 of Fig. 2 connected to an accelerometer 11 6, a pressure transducer 
1 1 8, and a delay circuit 120, all of which are included in the underwater seismic system 100. The 
accelerometer 1 16 and the pressure transducer 1 18 are physically located within the body 102 and 
respectively produce acceleration and pressure signal as described above. 
30 The underwater seismic system 1 00 includes a recorder 1 22 which is a conventional chart or 30 
drum recorder of a type having a repetitive sweep provided by sweep circuit 1 24. Each time the sweep 
circuit 1 24 initiates a sweep across the chart (not shown} a trigger signal is generated by a trigger 
signal generator 126. 

The delay circuit 1 20 receives the trigger signal, delays it in time, and then couples the trigger 
35 signal via the cable 1 04 to a sound source 1 28 located in the body 1 02. The delayed trigger signal 35 
activates the sound source 1 28 causing it to generate a sound pulse such as the pulse 1 08 of Fig, 3. 

Reflected pulse such as the pulse 11 2 of Fig. 3 are detected at intervals by hydrophones receivers 
1 30 which generate response signals indicating that reflected pulses have been detected. The 
amplitude and duration of each response signal will generally correspond to the amplitude and duration 
40 characteristic of each reflected sound pulse; however, variations in the characteristics of the successive 40 
response signals will normally not affect the operation of the underwater seismic system 1 00 in a 
material way. 

Each response signal is received by the recorder 1 22 which then causes a point to be plotted by 
an electrically activated pen (not shown) that regularly sweeps the recorder chart (not shown). In 

45 effect, the recorder mechanically computes and plots the time elapsed from the generation of the 45 
trigger signal to the receipt by the recorder 28 of the response signal. The portion of the charge 
recording so produced will also be indicative of the distance to the reflecting layer giving rise to the ? 
reflected sound pulse, since the average distance to the reflecting layer during the travel time of the 
sound pulse and reflected sound pulse will be proportional to the time elapsed, the constant of 

50 proportionality being the velocity of sound divided by two. In a similar manner, the recorder 122 50 
simultaneously plots from other response signals the distances to other reflecting layers. 

The extent of the delay between the initiation of the sweep by the sweep circuit 1 24 and firing of 
the sound source 128 depends on the vertical displacement signal produced by the filter 50. The delay 
circuit 1 20 receives the displacement signal and delays the trigger signal by a fixed period plus a 

55 varying period, the varying period being proportional to the displacement signal. Thus, as the towed 55 
body 1 02 is displaced downwardly, the delay in the firing of the sound source 1 28 is increased; as the 
towed body 102 moves upwardly, the delay is decreased and the firing is effectively advanced. 
Chanaes in the delav in firino of the sound source 128 will be related to changes in the vertical 
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velocity of sound in water. With such a constant of proportionality, the vertical displacement of the 
towed body 1 02 will appear to be the same at successive firings, for the purposes of the recorder 122. 

A more detailed description of apparatus comprising an underwater seismic system such as the 
underwater seismic system 100 may be obtained from the Hutchins patent referred to above. 

5 Claims 5 
1. A filter for use in determining the vertical displacement of a body travelling through water 
relative to a datum, the body having an accelerator meter with a predetermined accelerometer gain 
that produces a signal indicative of the vertical acceleration of the body, the body having a pressure 
transducer with a predetermined pressure transducer gain that produces a position signal, corrupted by 
1 0 pressure transducer noise, indicative of the vertical displacement of the body with respect to the 1 o 

datum, the filter comprising: 

(a) a first filter portion, connectable to the pressure transducer to produce a filtered position 
signal, the first filter portion having a transfer function which is substantially the product of a 
predetermined gain divided by the pressure transducer gain and a first real rational polynomial function 

15 oftheformP t (s)/Q{s); 15 

(b) a second filter portion, connectable to the accelerometer to produce a filtered acceleration 
signal, the second filter portion having a transfer function which is substantially the product of the 
predetermined gain divided by the accelerometer gain and a second real! rational polynomial function of 
the form sP 2 (s)/Q{s); and 

20 (c) a first summer connected to the first and second filter portions to produce a signal 20 
proportional to the sum of the first and second filter portion signals where: P,(s), P 2 (s) and Q{s) are 
predetermined real polynomial functions in the complex variable s; 

the degree of Q(s) in an integer n greater than 2, and the degree of each of P,(s) and P 2 (s) is less 
than n; 

25 the magnitude of P^sl/Qts) is less than 1 for s=jw when w is any radial frquency within the 25 
frequency range of the pressure transducer noise; and 

P,(s) and P 2 (s) satisfy the relationship P,(s)+S 3 P 2 (s)=Q(s). 

2. A filter as claimed in claim 1 in which the zeroes of Q(s) are such that a transfer function ?/Q(s) 
has a 20 n db/decade roll-off in a frequency range of the pressure transducer noise. 
30 3. A filter as claimed in claim 2 in which the zeroes of Q(s) are non-zero. 30 

4. A filter as claimed in claim 3, on which the zeroes of Q(s) lie in the closed left-hand complex s 
plane. 

5. A filter as claimed in claim 2 in which the zeroes of Q(s) lie in the open left-hand complex s 
plane. 

35 6. A filter as claimed in claim 2 in which the zeroes of Q(s) are negative, real numbers. 35 

7. A filter as claimed in claim 6 in which the degree n of Q(s) is 3. 

8. A filter as claimed in claim 7 in which the zeroes are substantially equal to one another. 

9. A filter as claimed in claim 8 in which: 



P^sJ^w^+Sw^s+Wo 3 ; 

40 M 40 
and 

Q(s)=(s+w 0 ) 3 ; 

where, w 0 is a predetermined positive real constant. 

1 0. A filter as claimed in claim 9 in which the first filter portion comprises: 

45 a first low pass filter, whose transfer function comprises two poles substantially at s=-w 0 , for 45 

producing a filtered position signal; 

a first band-pass filter, whose transfer function comprises a zero at s=0 and two poles 
substantially at s=— w 0 for producing a filtered position signal; 

a second low pass filter, whose transfer function comprises a pole substantially at s=— w^ 
50 connected to the first band-pass filter to produce a filtered, first band-pass filter signal; and 50 
a second summing means for producing a signal proportional to a weighted sum of the first and 
second low pass filter signals, and the first band-pass filter signal. 

11 . A filter as claimed in claim 1 0 in which the second filter portion comprises: 

a second band-pass filter whose transfer function comprises a zero at s=0 and two poles 
55 substantially at s=— w^ for producing a filtered acceleration signal; and 55 
a third low pass filter, whose transfer function comprises a pole at s=-w 0 , connected to the 
second band-pass filter to producing a filtered second band-pass filter signal. 

1 2. A filter for determining the vertical displacement of a body travelling through water, 
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constructed and arranged for use and operation substantially as described herein with reference to the 
accompanying drawings. 

1 3. A method of determining the vertical displacement of a body travelling through water, 
substantially as described herein with reference to the accompanying drawings. 
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